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PROBLEM  STATEMENT 


The  objective  of  this  research  project  is  to  provide  a  comprehensive  study  and  to  develop  new 
computational  methodologies  in  the  area  of  mechanics  and  control  of  constrained  deformable  bodies 
as  applied  to  large  scale  flexible  mechanical  systems.  The  main  objectives  of  the  project  can  be 
summarized  as  follows: 

(1)  Develop  a  classification  of  the  methodologies  used  in  the  dynami^and-eontrol  of 
flexible  multibody  systems  in  order  to  define  the  basic  assumptions  underlying  each 
method.  This  classification  will  be  used  to  establish  guidelines  for  selecting  the 
appropriate  formulation  for  a  given  application.  The  new  classification  will  be  based 
on  the  ability  and  accuracy  of  the  method  in  representing  the  rigid  body  dynamics, 
the  degree  of  accuracy  of  each  method  in  representing  the  coupling  between  the  rigid 
body  and  deformation  degrees  of  freedom,  and  the  basic  assumptions  underlying  each 
method. 

(2)  Use  the  new  classification  to  define  for  the  first  time  the  relationship  between 
different  methods  used  in  flexible  multibody  dynamics,  and  examine  the  effect  of  the 
linearization  resulting  from  the  use  of  the  incremental  finite  element  formulations  on 
the  definition  of  forces  in  flexible  multibody  dynamics. 

(3)  Develop  new  and  efficient  computational  finite  element  procedures  for  the  nonlinear 
dynamic  analysis  of  constrained  flexible  multibody  systems.  These  new 
computational  methods  will  lead  to  accurate  representation  of  the  large  displacements 
and  can  be  used  to  circumvent  the  limitations  of  existing  methods  used  in  the 
nonlinear  dynamic  simulation  of  constrained  flexible  multibody  systems. 

(4)  Examine  the  fundamental  relationship  between  the  deformable  body  coordinate 
system,  the  boundary  conditions,  the  resonance  conditions,  and  control  forces  in  the 
dynamics  of  constrained  deformable  bodies. 

(5)  Demonstrate  the  use  of  the  new  procedures  and  formulations  developed  in  this 
research  project  using  flexible  multibody  applications. 

The  objectives  of  this  research  project  have  not  been  modified  from  the  original  application. 


2.  MAIN  RESULTS 

Since  the  start  of  this  project  in  May  of  1997,  the  following  tasks  have  been  completed: 

1 .  A  review  of  past  and  recent  developments  of  the  dynamics  of  flexible  multibody 
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systems  is  presented  [1].  The  objective  was  to  review  the  basic  approaches  used  in 
the  computer  aided  kinematic  and  dynamic  analysis  of  flexible  mechanical  systems, 
and  to  identify  future  directions  in  this  research  area.  Among  the  formulations 
reviewed  are  the  floating  frame  of  reference  formulation,  the  finite  element 
incremental  methods,  large  rotation  vector  formulations,  the  finite  segment  methods, 
and  the  linear  theory  of  elasto-dynamics.  Linearization  of  the  flexible  multibody 
equations  that  results  from  the  use  of  the  incremental  finite  element  formulations  is 
demonstrated  and  discussed  [1,2]. 

2.  An  absolute  nodal  coordinate  formulation  is  developed  for  the  dynamic  analysis  of 
flexible  bodies  that  undergo  large  rotations  and  deformations  [1-8].  In  this  absolute 
nodal  coordinate  formulation,  no  infinitesimal  or  finite  rotations  are  used  as  nodal 
coordinates,  instead  global  displacements  and  slopes  are  used  as  the  nodal 
coordinates.  Using  this  set  of  absolute  coordinates,  exact  representation  of  the  rigid 
body  dynamics  can  be  obtained,  and  an  arbitrary  rigid  body  motion  leads  to  a  zero 
strain  in  the  finite  element.  The  absolute  nodal  coordinate  formulation  leads  to  a 
simple  constant  mass  matrix.  Furthermore,  using  this  formulation,  beam  and  plate 
elements  can  be  considered  as  isoparametric  elements. 

3.  The  absolute  nodal  coordinate  formulation  was  applied  to  several  flexible  multibody 
system  applications.  The  results  obtained  were  compared  with  the  results  obtained 
using  the  floating  frame  of  reference  formulation,  which  is  the  most  widely  used 
method  in  flexible  multibody  simulations.  Excellent  agreement  between  the  two 
methods  was  observed  in  the  case  of  small  deformation  problems  [9]. 

4.  The  effect  of  the  choice  of  the  element  coordinate  system  in  the  non-incremental 
finite  element  absolute  nodal  coordinate  formulation  [12-13]  was  investigated.  This 
formulation  leads  to  a  constant  and  symmetric  mass  matrix  and  highly  non-linear 
elastic  forces.  It  was  demonstrated  in  this  project  that  different  element  coordinate 
systems  which  are  used  for  the  convenience  of  describing  the  element  deformation 
lead  to  similar  results  as  the  element  size  is  reduced  [12].  In  order  to  demonstrate  this 
fact,  two  element  frames  are  used  in  this  study;  the  pinned  and  the  tangent  frames. 
The  pinned  frame  has  one  of  its  axes  passes  through  two  nodes  of  the  element,  while 
the  tangent  frame  is  rigidly  attached  to  one  of  the  end  points  of  the  element. 
Numerical  results  obtained  in  this  project  using  these  two  frames  are  found  to  be  in 
a  good  agreement  as  the  element  size  decreases. 

5.  The  equivalence  of  the  finite  element  formulations  used  in  flexible  multibody 
dynamics  has  been  the  subject  of  extensive  study  in  this  research  project.  This 
equivalence  was  used  to  address  several  fundamental  issues  related  to  the 
deformations,  flexible  body  coordinate  systems,  and  the  geometric  centrifugal 
stiffening  effect.  The  equivalence  of  the  inertia  forces  used  in  different  finite  element 
formulations  was  previously  established  [6].  The  relationship  between  the 
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coordinates  used  in  the  absolute  nodal  coordinate  formulation  and  the  coordinates 
used  in  the  floating  frame  of  reference  formulation  was  developed.  This  relationship 
is  used  to  obtain  the  highly  nonlinear  expression  of  the  strain  energy  used  in  the 
absolute  nodal  coordinate  formulation  from  the  simple  energy  expression  used  in  the 
floating  frame  of  reference  formulation.  It  was  shown  in  this  study  [12]  that  the 
source  of  the  nonlinearity  in  the  elastic  forces  is  due  to  the  finite  rotation  of  the 
element.  The  results  obtained  clearly  demonstrate  that  the  instability  observed  in  high 
speed  rotor  analytical  models  due  to  the  neglect  of  the  geometric  centrifugal 
stiffening  is  not  a  problem  inherent  to  a  particular  finite  elemenfTJtlfnulation.  Such 
a  problem  can  be  only  avoided  by  considering  the  known  linear  effect  of  the 
geometric  centrifugal  stiffening  or  by  using  a  nonlinear  elastic  model  as  recently 
demonstrated.  Fourier  analysis  of  the  solutions  obtained  was  used  to  shed  new  light 
on  the  fundamental  problem  of  the  choice  of  the  deformable  body  coordinate  system 
in  the  floating  frame  of  reference  formulation  [12, 13]. 

6.  Two  procedures  for  the  computer  implementation  of  the  non-incremental  absolute 
nodal  coordinate  formulation  were  proposed  for  flexible  multibody  applications  [10]. 
These  two  procedures  utilize  the  fact  that  the  mass  matrix  obtained  using  the  absolute 
nodal  coordinate  formulation  is  constant  in  order  to  obtain  an  optimum  sparse  matrix 
structure.  In  the  first  procedure,  a  QR  decomposition  of  a  modified  constant 
connectivity  Jacobian  matrix  is  obtained  for  the  deformable  body.  A  constant  velocity 
transformation  is  used  to  obtain  an  identity  generalized  inertia  matrix  associated  with 
the  second  derivatives  of  the  generalized  coordinates,  thereby  minimizing  the  number 
of  non-zero  entries  of  the  coefficient  matrix  that  appears  in  the  augmented 
Lagrangian  formulation  of  the  equations  of  motion  of  the  flexible  multi-body 
systems.  In  the  second  procedure,  a  Cholesky  decomposition  of  the  constant  mass 
matrix  is  used  to  define  another  constant  velocity  transformation  matrix.  This 
alternate  procedure,  which  has  the  same  computational  advantages  as  the  one  based 
on  the  QR  decomposition,  leads  to  a  square  velocity  transformation  matrix. 

7.  The  use  of  the  non-incremental  finite  element  absolute  nodal  coordinate  formulation 
in  the  analysis  of  large  deformation  and  rotation  problems  was  demonstrated  using 
several  structural  and  multi-body  applications  [8,  15].  The  results  obtained  are 
compared  with  the  results  obtained  using  existing  methods.  These  results 
demonstrate  that  the  absolute  nodal  coordinate  formulation  can  be  considered  as  the 
basis  for  developing  a  new  generation  of  general  purpose  flexible  multi-body 
computer  codes  that  can  be  used  in  the  analysis  of  small  and  large  deformation 
problems. 

8.  The  performance  of  the  new  absolute  nodal  coordinate  formulation  developed  in  this 
investigation  was  compared  to  the  incremental  finite  element  methods  using  several 
structural  and  mechanical  system  applications.  Numerical  results  showed  that  there 
is  a  good  agreement  between  the  results  obtained  using  the  new  method  and  the 


3 


results  obtained  using  the  widely  used  incremental  finite  element  methods  when 
structural  systems  are  analyzed.  On  the  other  hand,  several  large  rotation  and 
deformation  multibody  problems  showed  that  the  absolute  nodal  coordinate 
formulation  is  more  efficient  and  leads  to  a  better  accuracy  [14,15,19]. 

9.  New  and  simple  elastic  force  models  were  developed  using  the  absolute  nodal 
coordinate  formulation.  Four  new  models  that  employ  a  nonlinear  strain- 
displacement  relationship  were  presented.  These  new  force  modej^  were  used  to 
develop  a  general  purpose  flexible  multibody  computer  algorithm.  A  comparison 
with  the  model  previously  developed  shows  a  significant  reduction  in  computer  time 
and  improvement  in  the  accuracy  using  a  smaller  number  of  elements  [13,16,18]. 

1 0.  The  problem  of  the  geometric  centrifugal  stiffening  in  rotating  beams  was  examined 
in  details.  It  was  shown  that  the  absolute  nodal  coordinate  formulation  automatically 
accounts  for  the  geometric  centrifugal  stiffening.  This  is  an  important  feature  of  the 
new  formulation,  particularly  when  rotor  systems  are  analyzed.  It  was  demonstrated 
that  the  absolute  nodal  coordinate  formulation  can  efficiently  be  used  to  predict  the 
natural  frequencies  and  mode  shapes  of  rotating  beams.  Numerical  results  are 
compared  with  the  results  previously  published  [20]. 

11.  It  was  demonstrated  that  the  absolute  nodal  coordinate  formulation  can  lead  to  a 
constant  mass  matrix  for  three-dimensional  elements  by  increasing  the  number  of 
nodal  coordinates.  Using  this  important  property  of  the  new  formulation,  an  efficient 
procedure  for  the  solution  of  the  large  deformation  in  spatial  flexible  multibody 
systems  can  be  developed  using  Cholesky  coordinates,  as  previously  demonstrated 
in  the  planar  case  [17]. 
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